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Abstract 

f-^. ■ Spin polarized states in nuclear matter with an effective nucleon-nucleon inter- 

O ! action are studied for a wide range of isospin asymmetries and densities. Based on 

a Fermi liquid theory, it is shown that there are a few possible scenarios of spin 
ordered phase transitions: (a) nuclear matter undergoes at some critical density 
a phase transition to a spin polarized state with the oppositely directed spins of 
' neutrons and protons (Skyrme SLy4 and Gogny D1S interactions); (b) at some crit- 

ical density, a spin polarized state with the like-directed neutron and proton spins 
CN \ appears (Skyrme SkI5 interaction); (c) nuclear matter under increasing density, at 

first, undergoes a phase transition to the state with the opposite directions of neu- 
tron and proton spins, which goes over at larger density to the state with the same 
direction of nucleon spins (Skyrme SkI3 interaction). 
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The issue of spontaneous appearance of spin polarized states in nuclear matter is a 
topic of a great current interest due to its relevance in astrophysics. In particular, the 
scenarios of supernova explosion and cooling of neutron stars are essentially different, de- 
pending on whether nuclear matter is spin polarized or not. On the one hand, the models 



On 



O 



with the effective nucleon-nucleon (NN) interaction predict the occurrence of spin instabil- 
ity in nuclear matter at densities in the range from qq to 6^0 for different parametrizations 



of the NN potential [l]-[4] (qq = 0.16 fm _ ). On the other hand, for the models with the 
realistic NN interaction, the ferromagnetic phase transition seems to be suppressed up to 
densities well above £0 [5]-[7]. 

Here the issue of spin polarizability of nuclear matter is considered with the use of 
rS [ an effective NN interaction. The main objective is to study the possible scenarios of spin 

ordered phase transitions in dense nuclear matter with Skyrme and Gogny forces. In par- 
ticular, we choose Skyrme SLy4 effective interaction, constructed originally to reproduce 
the results of microscopic neutron matter calculations [8]. We utilize Skyrme SkI3 and 
SkI5 parametrizations as well, giving a correct description of isotope shifts in neutron-rich 
medium and heavy nuclei [9]. Besides, we employ Gogny D1S interaction, widely used in 
nuclear structure calculations. The basic formalism is presented in detail in Ref. [4]. We 
are interested in studying spin polarized states with like-directed and oppositely directed 
spins of neutrons and protons. One should solve the self-consistent equations for the co- 
efficients £00, £30, £03, £33 in the expansion of the single particle energy in Pauli matrices in 
spin and isospin spaces 

£oo(p) = £ o(p) + £oo(p) -/%), 6o(p) = £30 (p), (!) 

Cos(p) = £03(P) ~/i03, 6b(p) =£33(P)- 

Here £o(p) is the free single particle spectrum, /x 00 and /i 3 are half of a sum and half of 
a difference of neutron and proton chemical potentials, respectively, and e 00 , e 30 , e 03 , £33 
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are the Fermi liquid (FL) corrections to the free single particle spectrum, related to the 
normal FL amplitudes U (k), t/ 3 (k) by formulas 

eoo(p) = ^E^o(k)/oo(q), ^o(p) = ^X>i(k)/ 3 o(q), k = (2) 

eos(p) = ^E^(k)/ 03 (q), ess(p) = ^E^s(k)/3s(q). 

The distribution functions /oo> /o3) /30) /33) m turn, can be expressed in terms of the com- 
ponents £ of the single particle energy and satisfy the normalization conditions for the 
total density g n + g p = g, excess of neutrons over protons g n — g p = ag, ferromagnetic 
(FM) Qi — Qi = Ag^ and antiferromagnetic (AFM) (g n ^ + g P i) — (g n [ + g p f) = Ag-^ spin 
order parameters, respectively (a being the isospin asymmetry parameter, g^ = g n ^ + g p ^ 
and g[ = g n ^ + g p ±, with g n ^, g n ± and g p -[, g p ^ being the neutron and proton number den- 
sities with spin up and spin down). The quantities of interest are the neutron and proton 
spin polarization parameters IT„ = gwT ~ gwi , n p = gpT ~ gpi , characterizing spin ordering in 
neutron and proton subsystems. 

Fig. la shows the density dependence of the neutron and proton spin polarization 
parameters at zero temperature for SLy4 force. The main qualitative feature is that 
for SLy4 force there are only solutions corresponding to the oppositely directed spins of 
neutrons and protons in a spin polarized state. The reason is that for SLy4 force the FL 
amplitude Lq, determining spin-spin correlations, is repulsive for all relevant densities, 
while the FL amplitude U 3 , describing spin-isospin correlations, becomes quite attractive 
at high densities. The critical density of spin instability in symmetric nuclear matter 
(a — 0), corresponding to AFM spin ordering (Ag^ 7^ 0, Ag^ = 0), is g c ~ 0.33 fm -3 . It 
is less than the critical density of FM instability in neutron matter, g c 0.59 fm~ 3 . Even 
small admixture of protons to neutron matter leads to the appearance of long tails in the 
density profiles of the neutron spin polarization parameter near the transition point to 
a spin ordered state. As a consequence, a spin polarized state is formed much earlier in 
density than in pure neutron matter. 

As seen from Fig. lb, for SkI5 force, oppositely to SLy4 force, there are only solutions 
corresponding to the same direction of neutron and proton spins in a polarized state. In 
the case under consideration the FL amplitude U3 is repulsive for all relevant densities, 
while the FL amplitude U\ becomes quite attractive at high densities. For SkI5 force, a 
phase transition to the FM spin state in neutron matter takes place at the critical density 
g c R3 0.28 fm~ 3 . It is less than the critical density of spin instability in symmetric nuclear 
matter g c w 0.43 fm~ 3 , corresponding to FM spin ordering (Ag^ ^ 0, Ag^ = 0). There 
are no long tails in the density profiles of the neutron spin polarization parameter at large 
isospin asymmetry. In the given case, a small admixture of protons to neutron matter 
even leads to the increase of the critical density of spin instability. 

Fig. 2 shows the neutron and proton spin polarization parameters as functions of 
density at zero temperature for SkI3 force. There are two types of solutions of the self- 
consistent equations in symmetric nuclear matter, corresponding to FM and AFM ordering 
of neutron and proton spins. Due to proximity of FL amplitudes Lq and C/3, the respective 
critical densities are very close to each other (g c 0.910 fm~ 3 for FM ordering and 
g c ~ 0.917 fm~ 3 for AFM ordering) and larger than the critical density of spin instability 
in neutron matter (g c 0.37 fm~ 3 ). When some admixture of protons is added to neutron 
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Figure 1. Neutron and proton spin polarization parameters as functions of density at zero 
temperature for (a) SLy4 force and (b) SkI5 force. 

matter, the last critical density is shifted to larger densities and a spin polarized state with 
the oppositely directed spins of neutrons and protons appears. Under increasing density 
of nuclear matter, the neutron spin polarization continuously increases till all neutron 
spins will be aligned in the same direction. Protons, at first, become more polarized 
with density and their spin polarization is opposite to the spin polarization of neutrons. 
But, after reaching the maximum, spin polarization of protons decreases and at some 
critical density spins of protons change direction, so that the spin ordered phase with the 
like-directed spins of neutrons and protons is formed. Then, beyond the critical density, 
the spin polarization of protons is continuing to increase until the totally polarized state 
with parallel ordering of neutron and proton spins will be formed. Thus, for SkI3 force 
nuclear matter undergoes at some critical density a phase transition from the state with 
antiparallel ordering of neutron and proton spins to the state with parallel ordering of 
spins. With increasing isospin asymmetry, this critical density increases as well. Note that 
there are no long tails in the density profiles of the neutron spin polarization parameter 
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Figure 2. Same as in Fig. 1, but for SkI3 force. Also the curves, corresponding to FM and 
AFM ordering in symmetric nuclear matter, are shown. 

at large asymmetries. 

Now we present the results of the numerical solution of the self-consistent equations 
with the D1S Gogny effective force for symmetric nuclear matter. The main qualitative 
feature is that for the D1S force there are only solutions corresponding to AFM spin 
ordering and there are no solutions corresponding to FM spin ordering. In Fig. 3, it is 
shown the dependence of the AFM spin polarization parameter Ag^J g as a function of 
density at zero temperature. The AFM spin order parameter arises at density g « 3.8f?o 
for the D1S potential. A totally antiferromagnetically polarized state (Ag^Jg = 1) is 
formed at g « A.3go- The neutron and proton spin polarization parameters for the AFM 
spin ordered state are opposite in sign and equal to 

IT = -IT = ^ 

Q 

For comparison, we plot in Fig. 3 the density dependence of the AFM spin polarization 
parameter for the Skyrme effective forces SkM* and SGII being relevant for calculations 
at small isospin asymmetry. 

It is necessary to emphasize that different behavior at high densities of the interaction 
amplitudes, describing spin-spin and spin-isospin correlations, lays behind this divergence 
in calculations with different effective forces. These results clearly indicate the necessity 
to construct a new generation of the energy functionals with the properly constrained 
time-odd part at high densities. Probably, these constraints will be obtained from the 
data on the time decay of magnetic field of isolated neutron stars [10]. 
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Figure 3. AFM spin polarization parameter as a function of density at zero temperature for 
the D1S Gogny force and the SkM*, SGII Skyrme forces. 
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